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ABSTRACT 

Strong 7-ray emission from cocoons of young radio galaxies is newly predicted. Con- 
sidering the process of adiabatic injection of the shock dissipation energy and mass 
of the relativistic jet in active nuclei (AGNs) into the cocoon, while assuming ther- 
malizing electron plasma interactions, we find that the thermal electron temperature 
of the cocoon is typically predicted in ~MeV, which is determined only by the bulk 
Lorentz factor of the relativistic jet. Together with the time-dependent dynamics of 
the cocoon expansion, we find that young cocoons can yield thermal bremsstrahlung 
emissions at energies ~ MeV. 
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1 INTRODUCTION 



Relativistic jets in active galactic nuclei (AGNs) are widely 
believed to be the dissipation of kinetic energy of relativistic 
motion with a Lorentz factor of order ~ 10 produced at the 
vicinity of a supermassive black hole at the galactic center 
(Begelman, Blandford and Rees 1984 for reviews). The jet 
in powerful radio loud AGNs (i.e., FR II radio sources) is 
slowed down via strong terminal shocks which are identified 
as hot spots. The shocked plasma then expand sideways and 
envelope the whole jet system and this is so-called a cocoon 
(Begelman and CiofB 1989, hereafter BC89). The cocoon is 
a by-product of the interaction between AGN jets and sur- 
rounding intra-cluster medium (ICM). The internal energy 
of the shocked plasma continuously inflates this cocoon. Ini- 
tially, the existence of the cocoon is theoretically predicted 
by Scheuer (1974). 

The first clear evidence for an X-ray cavity was dis- 
covered in the center of the Perseus cluster of galaxies by 
Boehringer et al. (1993). The thermal ICM is displaced by 
the radio lobes which are composed of the remnants of the 
decelerated jet. Then the X-ray surface brightness in those 
regions are significantly deceased. These cavities correspond 
to the cocoons. Most of the X-ray cavities are associated 
with low power AGN jets (i.e., FR I radio sources). Re- 
cent X-ray observations of radio galaxies shows us a further 
evidences of these X-ray cavities (e.g., Fabian et al. 2000; 
Blanton et al. 2001). Another X-ray observational evidence 
of the cocoon is the non-thermal emission around radio lobes 
(e.g., Feigelson et al. 1995; Isobe et al. 2002; Croston et al. 
2005). In some cases, those non-thermal emissions are as- 
sociated with FR II radio sources. In any case, there is no 



direct evidence of thermal emissions coming from the dilute 
thermal plasma inside the cocoon. 

In this paper, we propose that "a cocoon of a young 
radio galaxy" as a new population of 7-ray emitters in the 
universe. Up to now, little attention has been paid to the 
evolution of thermal temperature and number density of the 
cocoon. Recently we have investigated the evolution of its 
temperature and number densities by taking the proper ac- 
count of mass and energy injections by the relativistic jet 
(Kino and Kawakatu 2005; Kawakatu and Kino 2006). We 
found that the cocoon remains constant temperature whilst 
the number density increases as a cocoon becomes younger. 
This leads to our new prediction of bright 7-ray emission 
from the young cocoon. 



2 COCOON INFLATION BY DISSIPATIVE 
RELATIVISTIC JET 

Here we consider the time-evolution of expanding cocoon 
inflated by the dissipation energy of the relativistic jet via 
terminal shocks. The adiabatic energy injection into the co- 
coon is assumed here. We will compare the source age t and 
a cooling time scale and check the consistency at the last 
part of §2. Note that the injection process of kinetic energy 
and mass into the cocoon is "continuous" during t. It is dif- 
ferent from the "impulsive" injection realized in gamma-ray 
bursts (GRBs) and supernovae. 

The time-averaged mass and energy injections from the 
jet into the cocoon, which govern the cocoon pressure Pc 
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and mass density pc are written as 



Pc{t)V,{t) 



2Jj(tMj(t), 



(1) 



(2) 



wliere 7c, 14, T^^ , Jj and Aj, are the adiabatic index of 
the plasma in the cocoon, the volume of the cocoon, the ki- 
netic energy and mass flux of the jet, and the cross-sectional 
area of the jet, respectively. The total kinetic energy and 
mass flux of the jet are T^^ = pjC^Fj^Uj, J, = pjTjVj where 
Pj, and Fj are mass density and bulk Lorentz factor of the 
jet (Blandford and Rees 1974). Hereafter we set uj = c. 
The total kinetic power of the relativistic jet is deflned as 
Lj = 2T^^ {t)Aj{t) and it is assumed to be constant in time. 
Although little attention has been payed to the mass in- 
jection Eq. ([21 up to now, it is of great signiflcance to take 
account of the Eq. for deriving the cocoon temperatures. 
Hence we take up the the Eq. ([2]) to evaluate the temper- 
atures of the cocoons. In contrast, the energy equation Eq. 
lU has been widely utilized in the literatures of the AGN 
bubbles in various ways (e.g., BC89; Dunn and Fabian 2004). 

The kinetic power dominance in the flow is postulated 
in this work in accordance with the observational indica- 
tions (e.g., Leahy and Gizani 2001; Isobe et al. 2002; Cros- 
ton et al. 2005). The jet is assumed to be cold since the hot 
plasma produced at the central engine usually cool down 
very quickly (e.g., Iwamoto and Takahara 2004). As for the 
mass and kinetic energy flux of powerful relativistic jets, nu- 
merical simulations tell us that no significant entrainment of 
the environmental matter takes place during the jet propa- 
gation (Scheck et al. 2002). According to this, the mass and 
kinetic energy flux of the jet are regarded as constant in 
time. Then, the conditions of T^^ = const., and Jj = const, 
leads to the important relations of 



pj{t)A^{t) — const, Fj(t) = const. 



(3) 



In fact, the constant Fj agrees with the relativistic hydro- 
dynamic simulations (e.g., Marti et al. 1997; Scheck et al. 
2002). In order to evaluate Lj, we use the shock jump con- 
dition of F?pj = /3hsPiCM (Kawakatu and Kino 2006) where 
Ais(= vhs/c) and picM is the advance speed of the hot spot 
/3hs = 10"^/3_2 (Liu et al. 1992; Scheuer 1995) and the mass 
density of ICM, respectively. Using, the jump condition, Lj 
is given by 

Lj = 2 X 10*"^ i?kpc/3-2n-2 erg s"' (4) 

where we use Aj{t) = iTR^^{t), and the hot spot radius 
-Rhs is given by -Rkpc = -Rhs(10^ yr)/l kpc. As a fidu- 
cial case, we set the number density of surrounding ICM 
as nicM(d) ~ picm(d)/mp — 10~^ cm~^n_2(rf/30 kpc)^^ 
where d is the distance from the center of ICM and cocoon. 
Since the change of the index from —2 does not change the 
essential physics discussed in this work, we focus on this case 
for simplicity. Since Lj is the ultimate source of the phenom- 
ena associated with the cocoon, all of the emission powers 
which will appear in §3 should be less than Lj. 

The number density of total electrons in the cocoon 
ne{t) is governed by the cocoon geometry and its plasma 
content. For convenience, we define the ratio of "the vol- 
ume swept by the unshocked relativistic jet" to "the vol- 



ume of the cocoon" as A(t). Hereafter we denote Vc{t) = 
2{-K/2,)Tl^Zl^{t), satisfies Zhs(t) = PhsCt, Rc, and 7^ = 
Rc/Zi^B < 1 as the cocoon volume, the distance from the 
central engine to the hot spot, is the radius of the cocoon 
body, and the aspect-ratio of the cocoon, respectively (e.g., 
BC89; Kino and Kawakatu 2005). Postulating that TZ and 
Zhs/Rhs are constant in time, A{t) is evaluated as 

2Ai{t)vit 



Ait) 



Vc{t) 



0.4 TZ ^ RkpcZ^o /3_l, 



(5) 



where Z30 = ^ha(10'^ yr)/30 kpc. Note that, in the case, 
the time dependence of A is deleted since K oc and 
Aj cc . We stress that this case satisfies the observational 
indication of ~ const (e.g., Conway 2002). Eq. (O tells us 
that A is of order unity. Actually it is seen in some numerical 
simulations (e.g., in Fig. 2 of Scheck et al. 2002). The cocoon 
mass density Pc{t) is controlled by the mass injection by the 
jet and it can be expressed as 

pc{t) ~ rjpj(t)^ 

= /3Lrj-VicM(^hs(t))^, (6) 

where we use the shock condition of F?pj = P'^^PiCM- Adopt- 
ing typical quantities of FR II sources (Begelman, Blandford 
and Rees 1984; Miley 1980; Bridle and Perley 1984), we ob- 
tain 

-2 



ne(t) ~ 4 X 10"=yin_2Fio/3: 



t 



IQT yr 



(7) 



where F = lOFio, and A = ^/0.4. Here we assume that the 
mass density of the pair plasma is heavier than that of 
electron-proton one, and then we adopt pc ~ 2mene in the 
light of previous works (Reynolds et al. 1996; Wardle et al. 
1998; Sikora and Madejski 2000; Kino and Takahara 2004). 
However, the mixture ratio of pair and electron-proton 
is still open. If we assume completely pure electron-proton 
content in the jet, too small is required and it conflict with 
that of non-thermal electrons (Kino and Takahara 2004). 

Let us estimate the electron (and positron) temperature 
(Te) and proton temperature (Tp). From Eqs. ([T]) and ((2)1 
together with the equation of state 

Pc « 2ncA;rc, (8) 
we can directly derive the temperatures as 



fcTc 



1 Fio MeV, kTp ^ 2 Fin GeV 



(9) 

where we adopt the two temperatures condition of kTe ~ 
{me/mp)kTp. It should be stressed that the temperatures 
are governed only by Fj . It is also worth noting that the geo- 
metrical factors in Eqs. ((l| and ([2} are completely cancelled 
out. Actually, the Fj dependence of Eq. ((9)l well coincide 
with the result of hydrodynamic simulations of relativistic 
outflows (Fig. 5 in Marti et al. 1997). One can naturally un- 
derstand these properties by comparing the well-established 
properties such as supernovae and GRBs. Constant tem- 
perature in AGN jet can be realized by the "continuous" 
energy injection into the expanding cocoon whilst tempera- 
tures of astrophysical explosive sources such as gamma-ray 
bursts and supernovae would be decreased because of "im- 
pulsive" injection of the energy. Since the shock dissipation 
of the relativistic flow into non-relativistic one, in general, 
requires the energy conversion of whole kinetic energy den- 
sity FjPjC'^ into internal one (Piran 1999). Thus the resultant 
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temperatures are uniquely governed by Fj and they remain 
to be constant in time. Similarly, in the studies of continous 
steller winds, the constant temperature has been predicted 
for a hot interior consist of the shocked wind (Weaver et al. 
1977). 

Here we examine the time scale of the Coulomb in- 
teraction between protons and electrons. The time-scale 
of energy transfer from the protons to electrons is given 
by tep ~ {npatc)~^ where at = 47r(e'^/fcre)^ In Ac is the 
transport cross section for electron-proton collision. The 
coulomb logarithm is written as In Ac ~ ln{3kTeXD/e^) 
where Ad = {kTe/innee^y^^ is the Debye length (Totani 
1998). A typical case of hot spots in AGN jets, we obtain 
In Ac ~ 50. Therefore, even using the maximal proton num- 
ber density rip ~ rie, tep satisfies 

^^5xl0^e?o.c-(^) (10) 

where Oe = kTe/nieC^ = 106 lo, and fie = ne(10'^ yr)/4 x 
10~^ cm~'^, are the electron temperature in unit of rUeC? , 
and normalized number density of thermal electrons, respec- 
tively. As mentioned before, recent studies suggest the ex- 
istence of large amount of pairs in AGN outfiows which 
lead to much smaller rip. Hence Eq. (|10|) shows the minimum 
value of tep{t)/t. Thus the energy transfer from protons to 
electrons is inefficient by the Coulomb coupling unless the 
cocoon is much younger than t ~ lO'' yr. 

Next we evaluate the time scale of thermal 
bremsstrahlung cooling in the cocoons. It is well known 
that thermal bremsstrahlung is inefficient for the dilute 
plasma since its emissivity shows oc nKt) where ne{t) 
is the electron number density in the cocoon. For the 
shock-heated electrons with the temperature of Qe ~ Fj, 
the cooling time of the bremsstrahlung per unit vol- 
ume is estimated as tbrom ~ FmeC^rie/ebrcm where the 
bremsstrahlung emissivity in the relativistic regime is 
ebrem = 1.3 X 10"^^ey^n^(l + 2.6Qe) erg cm"^ s"^ with 
the Gaunt factor of 1.2 (Rybicki and Lightman 1979). The 
condition of tbrcm(i) > t 

^^-5x10 e,o Ue (^^j (11) 

actually holds. Therefore most of the shock dissipation en- 
ergy is deposited into the cocoon without suffering strong 
radiative cooling and our treatment of adiabatic energy in- 
jection in Eq. ([1} is verified for t < ibrcm(t). We limit our 
attention on this case in the present work. 

On the thermalizations of electrons and protons, it is 
worth to refer the recent studies with the particle-in-cell 
(PIC) simulations, since we only examine the simple case 
of the classical Coulomb interaction. PIC simulations be- 
gin to shed light on the complicated microscopic dynam- 
ics with in the relativistic coUisionless shock. Using one- 
dimensional PIC simulations, Shimada and Hoshino (2000) 
revealed that the collision and merging processes among 
the coherent waves are accompanied by the strong thermal- 
ization of electrons. The results of three-dimensional PIC 
simulations (Nishikawa et al. 2003; Frederiksen et al. 2004) 
also tell us that electron populations are quickly thermal- 
ized whilst ion population tend to retain distinct bulk speed 
and thermalize slowly. In the result of Shimada and Hoshino 



(2000), (Fig. 4 in their paper), we see that the proton energy 
is transferred to the electrons, then the electrons are heated- 
up by protons. To sum up, PIC simulations imply the quick 
thermalization of electron populations by plasma waves and 
their associated instabilities such as the two-stream instabil- 
ity. Therefore, our estimate of Te would correspond to the 
lower limit of Te . It is not the purpose of this paper to derive 
more realistic Te in detail. 



3 EMISSIONS FROM A YOUNG COCOON 

3.1 Thermal MeV bremsstrahlung emission 

The time-dependence of the thermal bremsstrahlung lumi- 
nosity Lbrem iS giveU by I/brom(t) OC ne{t)Te ^^Vc{t) (X 

based on the cocoon expansion shown in the previous sec- 
tion. Hence it is clear that a younger cocoons are brighter 
bremsstrahlung emitters than older cocoons. In a similar 
way, brighter synchrotron luminosity has been expected for 
younger radio galaxies (Readhead et al. 1996; Begelman 
1996). With relativistic thermal bremsstrahlung emissivity 
(Rybicki and Lightman 1979), the luminosity of the opti- 
cally thin thermal bremsstrahlung emission uL^, at energies 
~ 1 MeV is estimated as 

Lbrom(t) ~ 2 X 10« n^Tl^eli'' (l^F^) (12) 

Here we omit the redshift (z) factor merely for simplicity. 
Eq. (|12p explains the reasons why for the no detection of the 
thermal emission from older cocoons. One is simply because 
it is not very bright. The other is because the predicted en- 
ergy range is ~ 1 MeV, the MeV-7 astronomy is still imma- 
ture and it is sometimes called as "sensitivity gap" compared 
with the energy range below 10 keV and above GeV ranges 
(Takahashi et al. 2004). 

For example, the bremsstrahlung emission from the co- 
coon located at a typical distance oi D = 10^ Mpc (O'Dea 
and Baum 1998) is examined here. In Fig. 1, we show the 
predicted values of uF,^ for the cocoons with t = 10^ yr 
and t = 10* yr. The cocoon with t = 10^ yr have uF,^ ~ 
lO"^'' erg cm'^^s"^. The detection threshold of SPI instru- 
ment on board the INTEGRAL satellite is about vFv 
10~^ erg cm~^s~^ at 1 MeV. For a young cocoon with 
t = 10* yr, the predicted luminosity is ^ 10^ times larger 
than that uFi, ~ 10~^^ erg cm~^s~'^. This is still less than 
the threshold of INTEGRAL. This may be the reason why 
no clear detection of MeV emission from young cocoons up 
to now. Fig. 1 shows that the XMM/Newton satellites can 
detect the low energy part of the thermal bremsstrahlung 
from young cocoons in principle. Hence some of extragalactic 
unidentified X-ray sources could be attributed as the low en- 
ergy tail of the bremsstrahlung emissions. Interestingly, the 
recent observation by XMM/Newton reveals that the spec- 
trum of young radio-loud AGN B1358-I-624 actually shows 
the power-law slope close to the bremsstrahlung's one (Vink 
at al. 2006). 

In MeV energy band, a proposed mission of detector 
SGD on board the NeXT satellite with the eye up to 
0.6 MeV (Takahashi et al. 2004) could detect the thermal 
MeV emission from those located slightly closer and/or 
younger with smaller Lorentz factor. Lastly, other future 
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Figure 1. Model prediction of MeV-peaked thermal 
bremsstrahlung emission from cocoons located at _D = 10^ 
Mpc. The predicted emission from young cocoon is brighter 
enough to detect in X-ray band whilst that from an old cocoon 
is much darker than the detection limits (Hasinger et al. 2001; 
Roques et al. 2003; Takahashi et al. 2004). 

mission, Advanced Compton Telescope (ACT), is worth 
to be noticed. The sensitivity of ACT is expected to 
be a significantly improved compared with INTEGRAL 
one (http:/ /heseweb. nrl.navy.mil/gamma/detector/3C 
/3C_sens.htm). Although they focus on the super- nova 
science (Milne et al. 2002) for the moment, the ACT would 
be a promising tool also for the young cocoon science. 

At a glance, one may think it hard to distinguish over- 
lapping emissions from the core of the AGN with limited 
spacial angular resolution of current satellites. Time vari- 
abilities of observed spectra is the key to distinguish them. 
It is obvious that the cocoon emission is constant in time 
whilst various emissions from the core of AGN should be 
highly variable. Hence steady emissions are convincingly 
originated in cocoons. Furthermore, the averaged spectral 
index of AGN core emissions at X-ray band (Koratkar 
and Blaes 1999; Kawaguchi et al. 2001) are softer than 
the bremsstrahlung emission discussed in the present work. 
Hence the difference of the spectral index is also a useful 
tool to figure out the origin of the emission. 

3.2 Non-thermal emissions 

Non-thermal emission from AGN jets is another key ingre- 
dient to investigate their physics. Non-thermal synchrotron 
emission from the radio lobes due to the relativistic electrons 
is well known characteristic of AGN jets (Miley 1980; Bridle 
and Perley 1984) . Recently, inverse-Compton (IC) emissions 
from large scale jets have been also intensively explored both 
theoretically and observationally (e.g., Celotti and Fabian 
2004; Croston et al. 2005). 

The properties of IC emissions from young cocoons is 
discussed here. For this purpose, we firstly consider the prop- 
erties of the synchrotron emission. The magnetic flux con- 
servation is assumed here during the jet propagation which 
is given by 

Bhs(t)-Rhs(t) = const. (1 ^ y ^ 2) 

where V is a parameter expressing the configuration of the 



magnetic field in the hot spot. The magnetic flux from the 
central engine is assumed to be constant in time. The case of 
constant Y — 1 shows the purely toroidal-dominated mag- 
netic fleld whilst F = 2 is relevant to the purely poloidal- 
dominated magnetic field. Using Y, the time dependence of 
the synchrotron luminosity at the hot spot I/hs.syn(i) may 
be given by Lhs.syn(t) oc RL(t)'y^Bl,{t)n'^'^ {^,1) oc t'^'^+i 
where we assume that the number density of the non- 
thermal electrons n^^(7, t) is proportional to n^{t), and 7 is 
constant in time because the synchrotron cooling time tend 
to be longer than the sound crossing time at the hot spot 
(e.g., KT04). Taking the observational fact of the large num- 
ber of the CSOs in spite of their young age, the larger syn- 
chrotron luminosity is required for younger sources (Begel- 
man 1996; Readhead at al. 1996). Qualitatively, the model 
well reproduce this observational properties of young radio 
galaxies. 

To evaluate IC emission of the cocoon, it may be use- 
ful to define the quantities of /ssc(t) = Uaac{t) /Usyn{t) and 
/ic/cmb(^) = t^ic/CMB(^)/f^syn(t), where Usyn, f/ssc, and 
t^ic/CMB, are the energy densities of synchrotron photons, 
those of synchrotron-self Compton (SSC), and those of IC 
scattering of the Cosmic-Microwave Background (CMB), re- 
spectively. Photons with larger density is the dominant seed 
photons for IC scattering. We denote that the IC lumi- 
nosities for synchrotron photons and the CMB as Lssc and 
iic/CMB, respectively. 

For Ussc{t) < UcMB, we see that /ic/CMB(i) oc 
which implies that younger cocoon produce less IC/CMB 
photons in contrast to the case of bremsstrahlung ones. It 
is of great importance to examine whether the predicted 
frequency of IC/CMB emission overlapping in MeV band 
or not. According to the standard diffusive shock acceler- 
ation, the acceleration time scale is estimated as tacc = 
(27r7mec5)/(e_Bhs) where ^ is the parameter characterizing 
the mean free path for the scattering (e.g., Drury 1983). 
The maximum Lorentz factor of electrons 7 can be ob- 
tained by the equating tacc to the synchrotron cooling time 
tsyn ~ (67rmeC^)/((jT7C_B^g). This shows the Bhs dependence 
of 7 as 

7(t)(xB-//^(t). (13) 

Assuming that the strength of magnetic field in the cocoon 
B{t) is proportional to the one in the hot spot, the maximum 
frequency of the IC/ICM emission i^ic/cmb 7^i^cMB can 
be estimated as 

J^IC/CMB(t) 1 X lO^M (^Y^) (1^) 

where we denote B{t) cx Bhs(i) oc R^^ (t) oc t^^ , 7(t) = 
10''74(t/10'' yr)^^^, and 74 = 7/IO''. The typical value of 7 is 
adopted from Blandford (1990). From this, one can find that 
J^ic/CMB(i) of young cocoons much smaller than 1 MeV = 
2 X 10^" Hz. 

In the case of Usyn{t) > L'^cmb, the behaviour of Usyn{t) 
is given by -Layn(t) oc cZl^{t)Usyn{t) oc esyn(i) Vc(f). From 
this we obtain L'syn(i) oc UB{t)t~^- The model predicts that 
fsBc{t) oc in time and younger cocoon yields more SSC 
photons. Using this, Lssc can be estimated as I/ssc(t) oc 
The maximum frequency of the SSC emission Vssc can be 
evaluated with Eq. ([T3|l . i^ssc ~ 7^i^syn ~ 1 x 10® 7*B Hz 
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can be evaluated as 

~ 1 X 10"-/tB-s (^Yo^) ^^^^ 

where we use Eq. ((13} and the typical value of B is set as 
B{t) = 10"^ B_5(t/10'^ yr)-'^ G based on Blandford (1990). 
Thus, it is found that non-thermal emissions from younger 
cocoon reside in much lower energy band than MeV. 



4 SUMMARY 

We have investigated the luminosity evolutions of AGN co- 
coons together with the dynamical evolution of expanding 
cocoon. Below we summarize the main results of the present 
work. 

• We newly predict the bremsstrahlung emission peaked 
at MeV-7 band as a result of standard shock dissipation 
of relativistic jets in AGNs. The temperatures of cocoon is 
governed only by the bulk Lorentz factor of the jet Fj. The 
electron temperature Te relevant to observed emissions is 
typically predicted in the range of MeV for Fj ~ 10. Con- 
stant temperatures of plasma in the cocoon can be realized 
because of the continuous energy injection by the jet with 
constant Fj. It should be emphasized that the constant be- 
haviour of AGN cocoon temperatures is different from the 
well known cases of gamma-ray bursts and supernovae. In 
these sources, the temperatures decrease in time because the 
energy injection time scales are much shorter than their ages. 
Since larger number densities of thermal electrons are pre- 
dicted for younger cocoons, brighter thermal bremsstrahlung 
emission than that of older cocoon is naturally expected. 

• Additionally, non-thermal IC emissions from young co- 
coons are also investigated. Importantly, in contrast to the 
case of MeV thermal emission, the typical frequency of SSC 
and IC/CMB emissions are predicted to be decreased for 
younger cocoon, since the maximum Lorentz factor of rel- 
ativistic electrons are decreased. Therefore the typical fre- 
quencies of IC from a younger cocoon are at much lower 
than MeV ranges. 
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